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Abstract 

The behavior of the adjoint Wilson line in finite-temperature, SU(2), lattice 
gauge theory is discussed. The expectation value of the line and the associated 
■ excess free energy reveal the response of the finite-temperature gauge field 

to the presence of an adjoint source. The value of the adjoint line at the 
critical point of the deconfining phase transition is highlighted. This is not 
calculable in weak or strong coupling. It receives contributions from all scales 
and is nonanalytic at the critical point. We determine the general form of the 
free energy. It includes a linearly divergent term that is perturbative in the 
bare coupling and a finite, nonperturbative piece. We use a simple flux tube 
model to estimate the value of the nonperturbative piece. This provides the 
normalization needed to estimate the behavior of the line as one moves along 
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the critical curve into the weak coupling region. 



Typeset using REVTgX 



I. INTRODUCTION 



The context of this paper is the finite-temperature, deconfining phase transition of pure- 
glue, SU(2), lattice gauge theory. Many of the properties of this transition are now well- 
established . It is second order and is in the universality class of the three-dimensional Ising 
model. The order parameter is the expectation value of the Wilson line in the fundamental 
representation. This is a measure of the response of the gauge system to the presence of a 
nondynamical, static, external quark source. The response of the system to the presence of 
sources in other representations has also been studied 0-§j] through the finite-temperature 
expectation value of Wilson lines in higher dimensional representations. The adjoint and 
3/2 have received most of the attention. 

The expectation value of the Wilson line at temperature T in representation J is the 
exponential of the corresponding excess free energy 



For the fundamental representation in the confined phase, the free energy is infinite, and 
the expectation value is zero. In the deconfined phase, the free energy is finite, and the 
expectation value is finite. It is believed that the infinite free energy in the confined phase 
is due to a tube of color electric flux with finite energy per unit length that extends from 
the source to infinity. 

Earlier work ]3],|| indicates that the adjoint line is nonzero at all finite temperatures. 
Furthermore, if t is the deviation of the temperature or of gauge coupling from the critical 
value, then the adjoint line has a singular behavior at the critical point that is determined 
by the specific heat exponent a |3|.[4| . 



In the confined phase, the adjoint source can be screened by a loop of flux of finite length. 
Nevertheless, the nonanalytic behavior of (Li) also reflects the diverging fluctuations of the 
flux at the critical point. 



(Lj) = e-W 



(1.1) 



(In) =a + b±\t 



(1.2) 
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From now on, we will be discussing the adjoint line and will drop the subscript "1". 

The purpose of this paper is to discuss the expectation value of the adjoint line in more 
detail. It focuses on the magnitude of (L) and especially upon its value at the critical point 
(L) c . The expectation value is a function of two variables Nt and g 2 or, equivalently, A and 
T. On the critical curve in either of these planes, (L) c is a function of a single variable which 
is essentially the cutoff. As the discussion will show, (L) c decreases rapidly as the cutoff is 
increased. This is an important limit on attempts to verify (|1.2|) closer to the continuum 
limit. It will be apparent that there is some interesting physics in producing the estimates 
for the magnitude of (L). One reason is that, as a matter of principle, neither weak nor 
strong coupling expansions can be used in a straightforward way to calculate (L) c . The 
critical point is just exactly the point at which the expansions must fail. Furthermore, the 
values for (L) given by the first terms in either expansion near the actual critical point are 
far from the values given by numerical simulation. For example, at Nt = 2 and 4/g 2 = 1.87, 
the Monte Carlo result is about 0.025 while first order strong coupling gives 0.009, and first 
order weak coupling gives —0.001. At Nt = 4 and 4/g 2 = 2.296, where the Monte Carlo 
result is about 0.01, first order weak coupling gives —0.63. 

To make progress, we will need a different approach. First, the general form of F will 
be established. With spatial cutoff a, coupling g at scale a, temperature T, and critical 
temperature T c , 

F = -c(g 2 )+T c d(T/T c ). (1.3) 
a 

The first contribution is the linearly divergent self energy. The function c has a power series 
expansion. Its first term is the 0(g 2 ) Coulomb self energy. In a free gauge theory, that is 
all there would be to F. Interactions give higher order corrections to the coefficient c. They 
also give rise to the finite function d, which contains effects from the nonperturbative scale 
T c ~ Aqcd- 

To estimate d, we will have to add some physical ideas and employ models that are not 
systematic approximations. In one way of looking at it, d includes the binding of a gluon to 
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the adjoint source to screen it. Alternatively, one can think in terms of the flux tube model 
IIU in which the phase transition is driven by the energy-entropy competition of the flux 
tube. In such a picture, let R be the return probability of the flux tube, and let A be the 
cost of attaching it to the line. Then we can write 

(L) c ~ AR. (1.4) 

In the strong coupling region, A and R can be estimated. 

Section II discusses the general form in ( |1.3|) , while Sec. Ill considers the weak coupling 
expansion of c(g 2 ). Section IV has discussion of the flux tube model and its surprisingly 
good predictions for the critical coupling. Finally, the main work is in Sec. V where the 
estimates of d(0) and d(l) are made. Section VI concludes the paper with a reiteration of 
the main points, which can also be stated now: Three methods with somewhat different 
approaches to the physics lead to equivalent results: 

(L) c = e -lNTC^)+d(l)] (L5) 

with 

c « (8/3) [0.351/(4//) + .294/ (4/g 2 ) 2 } (1.6) 

and 

e- dW ^0.1. (1.7) 



II. GENERAL FORM OF F 

We will concentrate our attention on the free energy F rather than on the line. The 
reasons for this will become apparent as the discussion proceeds. 
The first task is to establish the form 

F=-c(g 2 )+T c d(T/T c ). (2.1) 
a 



This can be written in a number of different ways depending upon whether one prefers a 
notation suggested by a lattice or continuum cutoff and whether one prefers to use only 
dimensionless variables. In dimensionless lattice units, the form is 

aF = c(g 2 ) + d(N Tc /N T )/N Tc . (2.2) 

The form in Q2.1|) is essentially the same as that derived by Dotsenko and Vergeles [0. 
They treated the case of a Wilson loop at zero temperature and showed that 

(W(C)) = e^ c) . (2.3) 

ip contains a term that is linearly divergent and proportional to the length of the loop C . 
The coefficient of the linear divergence is a finite function of the bare coupling g. It does 
not contain terms with logarithmic divergence. After this term, which is proportional to the 
length of the loop, is subtracted from ijj, the remainder is a finite function of the renormalized 
coupling g(fi). We assume that the proof that they have given at zero temperature applies 
at finite temperature also. Then the finite term in F has the form 

Td 1 (T,g 2 (fi),fi) = Td 2 (T/fi,g 2 (fi)). (2.4) 

When /i is chosen at T c , it becomes 

Td 2 (T/T c ,g 2 (T c )) = T c d(T/T c ), (2.5) 

since g 2 {T c ) is a fixed number of order one. We have chosen to extract a T c in the final form 
rather than a T because it seems likely that this contribution is not zero at T = 0, i.e. d(0) 
is finite. Our primary interest is in the critical value d(l). 

The linear divergence comes from the integration region where all of the points at which 
the gluons attach to the line come together, and the propagators have their T = short 
distance structure. Thus, the coefficient of the linear divergence is independent of T. If 
the renormalized charge is eliminated in favor of the scale T c , then c must be a function of 
the dimensionless combination aT c . However, this ratio of the cutoff and physical scales is 

6 



determined by the bare coupling g (the coupling at momentum scale 1/a). Thus, c is a finite 
function of the bare coupling as indicated in ( |2.1| ). 

For an indication of why F is the preferred quantity, recall that 

(L) = e' NTaF , (2.6) 

that c ~ g 2 , and that g 2 ~ 1/ ln(A^) on the weak coupling part of the critical line. Thus, it 
is more convenient to deal with F and then deduce the behavior of (L). 

At this point, we can see how our problem takes shape. Along the critical line, Nt c 
increases as we move into weak coupling. The function c has a perturbative expansion and 
contributes a decreasing factor to (|2.6|) . The function cI(Nt c /Nt) has contributions from 
short and long distance. However, on the critical curve, we need d(l), which is inherently 
nonperturbative. It contributes a cutoff-independent factor £(1) = e~ d ^ to (L) c . Thus, 
there are two factors. The first is calculable and gives a rapid decrease of (L) c with decreasing 
lattice spacing and increasing N T . The second factor is nonperturbative and independent of 
the cutoff. 

The first two terms in c can be estimated from the the calculations of Heller and Karsch 
||. To estimate d(l), we will make use of simple models and matching to the strong coupling 
region. 

III. WEAK COUPLING AND c(g 2 ) 

The calculations that Heller and Karsch || have done allow us to estimate the first 
two terms in c. We use the data in their Table 5. Although c is an essentially short 
distance quantity, it should be noted that there are some possible finite-volume problems. 
The finite- volume, constant mode of the gauge field cannot be handled straightforwardly in 
perturbation theory. It it has not been included in the calculations of Ref. ||. The volume 
dependence of the numbers that they have computed gives a rough estimate for the effect 
of the deleted mode and the finite-volume corrections. The spatial lattice sizes are 8 3 , 10 3 , 
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and 12 3 . The effect on the leading term in c is at the 2% level. However, it can be seen 
from their table that there is a substantial N$ variation in the 0{g i ) numbers. When they 
are used to compute the coefficient of g 4 in c, the finite-size effects persist. 

A closer analysis of the numbers at nearly fixed TL = N$/Nt shows that there are also 
significant finite-spacing effects. At infinite cutoff, the coefficient of the linear divergence is 
independent of TL. Thus, our best estimate of the 0(g 4 ) contribution to the coefficient will 
come from the largest Nt at fixed Ns/Nt but without regard for the particular, fixed value 
of that ratio. The 10 3 x 6 and 12 3 x 7 numbers combine to give the "0.294" in (|3.1| ) below. 
This could easily be off by 20%. However, an error of this size will not have a significant 
impact on the discussion to follow, in which other, larger uncertainties arise. This procedure 
gives our estimate 

c « (8/3) [0.351/(4//) + 0.294/(4//) 2 ]. (3.1) 

The factor 8/3 is needed to change from the fundamental to the adjoint representation 
of SU(2). In higher order, the relationship is not the simple multiplicative factor that it is 
here. The first term in c is the lowest order contribution to the linearly divergent Coulomb 
energy. It would be the whole answer if the theory were free. The interactions correct the 
coefficient of the linear divergence through interactions on the scale of the cutoff. This is 
the origin of the second term in ( |3.1| ). However, the interactions also correct the energy on 
all other scales. These are finite changes and are contained in the term d in F. Most of our 
attention will be focused on estimating these effects. 

IV. FLUX TUBE MODEL FOR THE PHASE TRANSITION 

Long ago Patel H proposed a flux tube model for the phase transition. In this picture, 
the phase transition is driven by the energy-entropy competition of the flux tube attached to 
a color source. This picture can be used to estimate the position of the phase transition. We 
find that it does surprisingly well. This suggests that the model is capturing some important 
physics and that it can be used for other estimates in later sections. 
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Let the length of a basic unit of the flux tube be a. In the strong coupling region, a is 
the lattice spacing a. In weak coupling, it is of the order of the diameter of the flux tube, 
which in turn is of the same order as the nonperturbative scale 1/Aqcd and the inverse root 
string tension l/y^cr). We will use a — l/y/(a) in the weak coupling region. 

Each unit of flux, has an associated energy aa and Boltzmann weig ht e~ ad/T . The 
entropy factor for each unit is 7. For a free random walk of flux in three dimensions, there 
are six possible directions for each step, so 7 = 6. However, the flux tube has repulsive self 
interactions at crossings that suggest the use of a self-avoiding walk. In that case, a simple 
guess for 7 is 5. Numerical simulations of self-avoiding walks give 7 ~ 4.68 f9|. We will use 
that value. In other work ||, we have shown that the fractal dimension of the flux tube at 
T c is near the dimension of a self-avoiding walk. 

The phase transition occurs when the temperature has increased to the point that the 
product of the energy and entropy factors for each step is equal to one. 

In strong coupling, where aa 2 « Ing 2 , that gives 

(9 2 c ) Nt = 7 (4.2) 

or 

4/2 c 2 = 4(l/ 7 ) 1/7Vr . (4.3) 



For N T = 2, this gives 4/g 2 = 1.85. The Monte Carlo value is 1.87 |IJ. For N T = 1, QO| 
and Monte Carlo give 0.855 and 0.873 |pT| , respectively. 
In weak coupling a rearrangement of ([4.2|) gives, 



In 7 = aa/T c = ^(a)/T c . (4.4) 

A Monte Carlo calculation []T2[ of the ratio on the right gave 1.45. On the other hand, the far 
left term is 1.54. This is acceptable agreement given the rough nature of the estimates. For 
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Nt = 3, the phase transition is at 4/g 2 = 2.177 [I2[, which is just into the weak coupling 



region. At this coupling, the square root of the string tension is 0.492 [|T^]. This gives 
Nt c \/{(J(x 2 ) = 1.48, which also compares reasonably with In 7 = 1.54 in ( [4. 4| ) . 

The crossover between the strong and weak coupling regions is around g 2 ~ 4/g 2 w 2. 
Thus, the = 2 phase transition is in the strong coupling region but near the crossover. 
At 4/g 2 = 1.87, aa 2 = Ing 2 gives l/^/(aa 2 ) ~ 1.15 so that the lattice spacing is only slightly 
larger than the length l/^(a), which is the diameter of the flux tube. Thus, at this point, 
the lattice cutoff is not eliminating a great many flux tube configurations. A reasonable 
match of strong and weak coupling quantities in this region seems possible. 

One may worry that the strong coupling flux tube picture of the phase transition that 
uses the lowest order result aa 2 = \ng 2 for the T = string tension is too rough. However, 
the first correction to the string tension comes from graphs with four extra plaquettes. This 
is of relative order (g 2 )~ 4 . Even at the crossover g 2 »s 2, this is just 1/16. 

At this point, one may fairly inquire into the physical effects that have been neglected in 
the flux tube picture. The model asserts that the important finite temperature fluctuations 
of the gauge field are large scale fluctuations in the configuration of the flux tube that is 
attached to an external source. The flux tube is treated as a structureless object that has 
no internal excitations when the temperature is increased. As indicated, this may be a 
good approximation in strong coupling. In weak coupling, where the flux tube has a finite, 
nonperturbative size, this is a simplification that could be incorrect. In addition, the picture 
does not explicitly include the thermal fluctuations that are closed loops of flux or the more 
mundane excitations of short wavelength modes. The latter may not be too important since 
they are not highly excited at T = T c . The loops of flux that do not intersect the flux tube 
connected to the line cancel in the ratio that is the excess free energy of the static source. 
Part of the effect of those that do intersect may be already included in using the self-avoiding 
walk value 7 = 4.68. 
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V. ESTIMATES OF d(0) AND d(l) 



A. d(0) and the glueball mass 

One rough estimate of d(0) comes from the glueball mass. Let us view the glueball as 
a bound state of two gluons, each of which carries the adjoint representation of the gauge 
group. Then the glueball mass is an indication of the energy associated with the response of 
the vacuum to an adjoint source. The problem, of course, is that the gluons are "massless" 
while the line is infinitely massive. For nonrelativistic bound states, the system with one 
infinitely massive constituent has twice the binding energy and one half the radius of the 
equal mass system. On the other hand, a simple bag model calculation that includes only 
volume and kinetic energies gives a glueball mass a factor of 2 3 / 4 larger than the system with 



one infinitely massive constituent |I3|- If the effects of vacuum fluctuations and attractive 
color hyperfine interactions were included, they would act oppositely so that the ratio would 
not change greatly. This suggests that d(0) might be within a factor of two of the ratio 



of the glueball mass to the critical temperature. The latter has been measured at 5.6 [12 
Thus, for a rough guess with large uncertainty, put 

d(0) « 5.6. (5.1) 



B. Strong coupling estimate of d(0) 

It is easy enough to get the lowest order strong coupling contribution to aF(0). It comes 
from the strong coupling diagram that is a 1 x 1 x 1 tube of plaquettes in the inverse 
temperature direction: 

oF(0) = 41n# 2 . (5.2) 

In units of T c , this is 

F(0)/T c = 4iV Tc ln2 2 . (5.3) 
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On the other hand, we know that 



(g 2 ) NTc = 7 (5-4) 



so 



F(0)/T c = 4hi7^6.2. (5.5) 
At this level of approximation, the result is independent of g 2 . Notice that it is similar to 



the number in (|5.1| ). In the strong coupling region, there is nothing that corresponds to the 



short distance effects in c. Thus, the strong coupling F(0)/T c can also be thought of as d(0). 

C. d(l) from simple matching 

Recall that we are looking for an estimate of F(T C ) that is useful for decreasing g 2 and 
increasing Nt- This means that we need the value of d(l) in aF(T c ) = c(g 2 ) + <i(l)/iVy c . 

1. Monte Carlo 

The simplest thing to do is to determine F(T C ) from a Monte Carlo measurement of 
(L) c at one point in the weak coupling region where c can be reliably computed. With F 
and c given, d is determined. Unfortunately, there are no numerical experiments that have 
carefully looked at the critical value of (L) with small coupling and large Nt- It is easy to 
see from the published Monte Carlo data or from our discussion that (L) c will be small and 
somewhat difficult to measure at large Nt- At both Nt = 1 and 2, (L) c ps 0.025 |||14). The 
values at Nt = 4 and 6 are smaller and more uncertain — around 0.01 at Nt = 4 |15| and 
0.001 to 0.0014 at Nt — 6 Q. If, for the sake of discussion, we use these Nt = 4 and 6 
numbers, then the corresponding values for d(l) are 2.38 and 3.79 to 3.45. The values for 
£(1) are 0.093 and 0.022 to 0.031. 
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2. Strong coupling form 



The next simplest thing to do is to match the weak coupling form ( j2.2| ) to a strong 



coupling form at the crossover around g 2 = 2. As indicated in the introduction, it is not 
possible to calculate at T c in a straightforward expansion. Instead, we will employ a very 
simple model. It is inspired by the flux tube picture of the transition. We assume that the 
most important contribution at T c is captured in a very limited subset of the strong coupling 
diagrams. These diagrams are certain generalized cylinders of plaquettes in the fundamental 
representation. They are invariant under translations along the inverse temperature axis. 
On each three-dimensional slice, they consist of a self- avoiding loop that begins and ends at 
the location of the adjoint source. This model goes a bit too far in that the flux tube has 
repulsive self interactions but is not strictly self- avoiding. However, in a related context, we 
have found || that the fractal dimension of any field theory walk with repulsive interactions 
and interactions with background loops is close to the dimension 1.69 of a strictly self- 
avoiding walk. So perhaps some of the effects that have been neglected are partly accounted 
for in the assumption of a self-avoiding configuration of flux. 

With this restriction on the allowed graphs, it is possible to sum up their contribution 
at T c . We did this in the following way: The asymptotic form for the number of closed, 
self-avoiding configurations is [§] 



In addition, each configuration is weighted by the cost of a strip of N T plaquettes raised to 
a power that is the length of the path 



N(i) ~ r 7 /y . 



(5.6) 




(5.7) 



It follows that the critical point is at 




(5.8) 



as indicated previously. 



13 



The sum to be done is 

oo 

R = T,m)(ihY- (5-9) 

1=4 

For the terms through I = 12, we used the explicit numbers for N(l) given in JllJ and, for the 
rest, the asymptotic form (|5.6|) . The multiplicative constant implicit in Q5.6Q was adjusted 
to match the known value at I = 12. This gives R = 0.078. To get the value for (L) c , we 
must divide this by the "3" that relates the Wilson line to the adjoint character. The result 
is 

(L) c « R/3 w 0.026. (5.10) 
Notice that this agrees with the Nt = 1 and 2 Monte Carlo results. 

3. Matching 

Now select d(l) so that as g 2 approaches the crossover region from the weak coupling 
side, then (L) c approaches 0.026. The equations ( |2.2|) and ( |3.1|) are used. Using Nt = 2, 
A/g 2 = 1.87, and (L) c = 0.026, gives d(l) = 2.20. 

The sensitivity of d(l) to the details of the matching will be discussed in Sec. V.F. For 
now, we will accept 

d(l)«2.2 (5.11) 

as a reasonable estimate with considerable uncertainty. This can be compared with the 
estimate d(0) = 5.6 that came from the glueball mass. As one would guess, it appears that 
d(0) and d(l) are of the same order of magnitude. 

D. Return probability approach to d(l) 

This is an alternative approach to the flux tube picture of the phase transition. It 
eventually leads to the same calculation as did the previous method. There is a loop of flux 
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attached to the adjoint source. As the critical point is approached, this flux has increasing 
fluctuations. The number 0.026 computed above can be thought of as one third of the sum 
over lengths of the return probabilities of self-avoiding configurations of the flux tube. 

In Sec. V.C.2, we discussed this in the strong coupling region. To use the same picture 
in weak coupling, an elaboration is required. There is a cost associated with attaching the 
flux tube to the point source. The diameter of the flux tube is of the order of the inverse 
root of the string tension, which is many lattice spacings in the weak coupling region. The 
source is a point on the lattice. This mismatch of scales has an associated cost that we will 
call A. With this in mind, we consider the form 

(L) = ~AR. (5.12) 

Again, the 1/3 is the conventional factor in the relation between the line and the adjoint 
character. 

In the strong coupling region, the diameter of the flux tube is smaller than a lattice 
spacing so that A = 1, explicitly, in the lowest order strong coupling model described above. 
The first correction to A = 1 is 0{Nx{g 2 )~ & ) (when the line excludes two touching boxes on 
the sheet), which is small for the strong coupling values of Nt and g 2 that are relevant. A 
is unknown in weak coupling because it contains contributions from scales from the lattice 
spacing on up to the the size of the flux tube. However, as the crossover is approached these 
scales converge and A — > 1. 

Now consider the return probability R in the weak coupling region at the phase transition. 
Loops of arbitrarily large size contribute to R. This is large distance physics. We assume that 
strong coupling gives an adequate description of large distances. Thus, we use R = 0.078 in 
both regions. 

To determine d(l), use Eqs. ( |2.2j ), (|2.6j), and ( pTl| ) on the left hand side of ( |5.12j ) and 



evaluate it in the crossover region with A ~ 1. Evidently, this leads to the previous result. 
However, the discussion has provided a bit more insight into the nature of the physical 
approximations that are involved. 
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E. d(l) from matching F(0) 



This is another way to look at the physics of the problem. We assume that for all 
couplings, the change in the free energy as the temperature is increased from zero to T c is 
mainly due to the increasing large scale fluctuations of the flux tube and further that these 
large scales are adequately described by our model of the strong coupling flux tube. This 
will give us d(l) — d(0). To get d(0), we require continuity of the zero temperature free 
energy F(0) through the crossover region. Although this approach is significantly different 
from the two already discussed, it will turn out to be nearly equivalent in its net result. 

In the first step, we adjust d(0) so that the weak coupling expression for F(0) approaches 
the strong coupling value as g 2 is increased. Specifically, we require 

c/a + T c d(0) = (if a) In g 2 (5.13) 

or 

d(0) = 4N To lng 2 -N Tc c(g 2 ) (5.14) 

at the crossover. For g 2 = 2 and N T = 2.3 |7|, this gives d(0) = 4.85. It can be compared 



directly with the value d(0) rs 5.6 obtained from the glueball mass. Given the uncertainty 
in the latter, these are not inconsistent. 

The second step asserts that the change in the free energy 

[F(T C ) - F{0)]/T c = d(l) - d(0) (5.15) 

in weak coupling can be approximated by the strong coupling flux tube value. Our model 
for the flux tube gives a value for the change in the free energy that becomes independent 
of g 2 on the strong coupling side of the crossover. 

[F(T C ) - F(0)]/T c = - ln[(L) c ] sc - 4iV Tc \ng 2 (5.16) 

= -ln[(L) c ] sc -41n 7 (5.17) 

= -ln(0.026) -41n(4.68) (5.18) 

= -2.52 (5.19) 
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and 



= d(0) + [d(l) - d(0)] « 4.85 - 2.52 « 2.3. (5.20) 

This is not quite the same as the value 2.2 in ( |5.11| ). The method of the present section 
would give precisely the result of (|5.11|) if the values for ANt c In g 2 in ( |5.14j ) and (|5.19|) had 



been the same. In ( |5.14| ), we worked at the crossover for the best estimate of d(0). However, 



in the spirit of the method, a value well into the strong coupling region should be used 



in ( |5.19| ). So as a matter of principle, this method is not the same as that of Sec. V.C. 
However, 4Nt c In g 2 is a slowly varying function over the range from the crossover and into 
strong coupling. Thus, in practice the values for d(l) end up rather close, The difference is 
not significant considering the rough nature of our approximations. 

F. Sensitivity to the matching point 



There is some freedom in the details of the matching in the crossover region. We will 
consider these for the procedure of Sec. V.C. 3. Similar considerations apply to the other 
methods. Since all of the methods are equivalent or nearly equivalent to that of Sec. V.C. 3, 
it is sufficient to discuss that case. 

In Table |, we have chosen representative values for Nt, 4/g 2 , and (L) c . For the noninte- 
ger values of N T , 4/g 2 was calculated from ((Of). The results for d(l) and £(1) follow from 
( |2.2|) , ( |2.6|) , and (|3.1| ). We conclude that £(1) is somewhat sensitive to the details of the 



matching. However, except for the Nx = 3 cases, for which the (L) c values were arbitrarily 
chosen, the results for d(l) and C(l) are around 2.1 to 2.5 and 0.09 to 0.12, respectively. 

VI. NET RESULT 

While each of the methods end with essentially the same result, the physics is viewed 
from a different perspective in each of the cases. This helps us to understand the true nature 
of the approximation. 
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The final expression is 



IL) = e -[NTc(^)+(T c /T)d(T/T c )] ^ 

= e- NT < 9 ^C(T/T c ). (6.2) 



The function c is given in ( |3.1| ). At the critical point, 

(L) e = e- NTc{92) C(l) (6.3) 

with 

£(1) » 0.09 to 0.12. (6.4) 

This is several times higher than the Nt = 6 Monte Carlo numbers that were mentioned 
in Sec. V.C.I. However, we feel that those numbers must be considered tentative until 
confirmed by a simulation that is specifically focused on (L) c . The main technical weakness 
of our work is the expression for c(g 2 ). To reduce the uncertainties, it would be useful to 
have an infinite volume, 0(g 6 ) calculation of c. 

ACKNOWLEDGMENTS 

P.V. would like to thank U. Heller for useful discussions. This research was supported 
by the United States Department of Energy. For P.V., the grant numbers are DE-FG05- 
85ER250000 and DE-FG05-92ER40724. 



18 



REFERENCES 



[1] The seminal theoretical work is in B. Svetitsky and L. Yaffe, Nucl. Phys. B210[FS6], 
423 (1982); B. Svetitsky Phys. Rep. 132, 1 (1986). There is now a large literature of 
Monte Carlo results. Some recent references are: J. Engels, J. Fingberg, and M. Weber, 
Nucl. Phys. B 332, 737 (1990); J. Engels, J. Fingberg, and D. Miller, Nucl. Phys. B 
387, 501 (1992); J. Engels and V. Mitrjushkin, Phys. Lett. B 282, 415 (1992). 

[2] P. Damgaard, Phys. Lett. B 194, 107 (1987); K. Redlich and H. Satz, Phys. Lett. B 
213, 191 (1988). 

[3] J. Kiskis, Phys. Rev. D 41, 3204 (1990). 

[4] J. Fingberg, D. Miller, K. Redlich, J. Seixas, and M. Weber, Phys. Lett. B 248, 347 
(1990). 

[5] A. Patel, Nucl. Phys. B 243, 411 (1984); F. Takagi, Phys. Rev. D 34, 1646 (1986). 

[6] J. Kiskis, R. Narayanan, and P. Vranas, in The Rice Meeting, DPF90 edited by B. 
Bonner and H. Miettinen, (World Scientific, Singapore, 1990); J. Kiskis, R. Narayanan, 
and P. Vranas, in The Fermilab Meeting, DPF92 edited by C. Albright, P. Kasper, 
R. Raja, and J. Yoh, (World Scientific, Singapore, 1993); J. Kiskis, R. Narayanan, 
and P. Vranas, University of California at Davis, Report No. UCD-93-06, 1993 (to be 
published) . 

[7] V. Dotsenko and S. Vergeles, Nucl. Phys. B 169, 527 (1980). 
[8] U. Heller and F. Karsch, Nucl. Phys. B 251, 254 (1985). 

[9] M. Barber and B. Ninham, Random and Restricted Walks (Gordon and Breach, New 
York, 1970). 

[10] M. Okawa, Phys. Rev. Lett. 60, 1805 (1988). 

[11] R. Ben-Av, H. Evertz, M.Marcu, and S. Solomon, Phys. Rev. D 44, 2953 (1991). 

19 



[12] J. Fingberg, U. Heller, and F. Karsch, Nucl. Phys. B 392, 493 (1993). 

[13] This follows from the simplest version of the model developed in A. Chodos, R. Jaffe, 
K. Johnson, C. Thorn, and V. Weisskopf, Phys. Rev. D 9, 3471 (1974); A. Chodos, R. 
Jaffe, K. Johnson, and C. Thorn, Phys. Rev. D 10, 2599 (1974); R. Jaffe and J. Kiskis, 
Phys. Rev. D 13, 1355 (1976); T. DeGrand, R. Jaffe, K. Johnson, and J. Kiskis, Phys. 
Rev. D 12, 2060 (1975); R. Jaffe and K. Johnson, Phys. Lett. B 60, 201 (1976). 

[14] P. Vranas, Nucl. Phys. B(Proc.SuppL) 20, 337 (1991). 

[15] K. Redlich and H. Satz, Phys. Lett. B 213, 191 (1988). 

[16] M. Fisher and M. Sykes, Phys. Rev. 114, 45 (1959). 

[17] The value Nt = 2.3 for g 2 = 2 was obtained by interpolating between the Monte Carlo 
results at integer values of Nt- Eq. fl4.2| ) gives essentially the same result. 



20 



TABLES 



TABLE I. Matching points. 



N T 



(Li 



d(l) 



C(l) 



2.0 
2.2 
2.3 
2.4 
3.0 
2.0 
2.2 
2.3 
2.4 
3.0 
3.0 



1.87 

1.98 

2.0 

2.10 

2.177 

1.87 

1.98 

2.0 

2.10 

2.177 

2.177 



0.026 
0.026 
0.026 
0.026 
0.026 
0.020 
0.020 
0.020 
0.020 
0.020 
0.010 



2.20 
2.17 
2.12 
2.15 
1.86 
2.46 
2.43 
2.38 
2.42 
2.13 
2.82 



0.11 

0.11 

0.12 

0.12 

0.16 

0.085 

0.088 

0.092 

0.089 

0.12 

0.060 
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